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ABSTRACT 



In this work, we study the timing and spectral evolution of the transient low mass X- 
ray binary IGR J17091-3624 during its 2011 outburst. We present results obtained from 
observations with two instruments, Rossi X-ray Timing Explorer (RXTE) proportional 
counter array (PCA) and SWIFT X-ray telescope (XRT), between 03 February 2011 
and 10 April 2011. Apart from the detection of increase in X-ray flux while moving 
from the canonical low hard state to the hard intermediate state, we report the details 
of a transition from a low irregular variability state, commonly seen from most of the 
black hole X-ray binaries (BHXBs), to a regular, repetitive and large variability state, 
seen only from GRS 1915+105 ('p' class variability). By comparing the timing and 
spectral properties of IGR J17091— 3624 in the variable state with those of the '/)' class 
of GRS 1915+105, we conclude that till now the former is the only known source, which 
somewhat behaves like GRS 1915+105, and could be used as a missing link between GRS 
1915+105 and other canonical BHXBs in understanding their accretion and radiation 
mechanisms. 

Subject headings: accretion, accretion disks — black hole physics — X-rays: observa- 
tions — X-rays: binaries — X-rays: individual: (IGR J17091-3624, GRS 1915+105) 



1. Introduction 



IGR J17091— 3624 was detected as a bright variable transient source, with flaring activity 
in October 1994 (Mir/KVANT/TTM ), September 1996(BeppoSAX/ WFC), Septemb e r 200 1 (Bep- 



poS AX / WFC: (lln't Zand et al 



20031 )^. April 2003fINTEGEAL/IBIS: iKuulkers et al.l (|2003l)^. and 



July 2007 (|Capitanio et al.ll2009l l. In April 2003, IBIS, the gamma-ray imager on board the INTE- 



GRAL satellite detected the transient X-ray binary IGR J17091-3624 at the position of 17^09. r 



(R.A . (J2000.0)), -36°24.63' (declination (J2000.0)) while monitoring the galactic centre (|Kuulkers et al 



20031 ). A transition between soft and hard states and the presence of blackbody component 



during source softeni ng were observed durin g a joint spectral analysis of INTEGRAL /ISIE and 
RXTE/FCA data bv ICapitanio et al.l (|2006l ). They found that the spectral be havior and time 



variability of the source was similar to the black hole candidate 111743—322. iLutovinov et al 



(j2003l ) found that the power density spectrum (PDS) of IGR J17091-3624 can be well described by 
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the standard model of band-limited noise with the break frequency of 0.31 it 0.04 Hz and fitted 
the energy spectrum with the powerlaw model having powerlaw in dex of 1.43 zfc 0.03, commonl y 
observed in X-ray binary systems in the low hard spectral state (jSunvaev &: RevnivtsevI |2000| ) . 
SWIFT/XRT observations of IGR J 1709 1-3624 during the outburst in 2007 revealed another state 
transition from hard to soft, and refined the X-ray position of this source, which was consistent w ith 
the position of a compact radio counterpart in the VLA archived data (jCapitanio et al.ll2009l ). as 
well as the position of the infrared counterpart observed with the ESO New Technology Telescope 
(NTT) (jChatv et al.ll2008l ). Radio emission from this source with steep synchrotron e mitting spec- 



trum had been confirmed with the observations from the Very Large Array ( VLA) (jRupen et al 



20031 ) and the Giant Meterwave Radio Telescope (GMRT) (jPandev et al.ll2006l ). In spite of several 



efforts, little is known about the compact nature of this source. There is no optical observation of 
its binary counterpart. The renewed activity of the source was first detected by the SWIFT/HAT 
when the 15.0-50.0 keV hard X-ra y flux increased upt o 60 mCrab on 03 February 2011, compared 



its low hard level of ~ 10 mCrab (jKrimm et al.ll201ll ). This i s next to the large s t outb urst (~ 70 



mCrab) in 2007. Later using the SWIFT/X.BT observation, iKrimm &: Kenneal (|201ll ) confirmed 
that the outburst occurred in the source IGR J 1709 1-3624, and not in IGR J 17098-3628, which 
is also an X-ray transient 9.6' away from IGR J17 091-3624. The INTEG RAL/IBIS also detected 



the source simultaneously with the SWIFT/XJIT (jPel Santo et al.ll201ll ). During 2011 outburst. 



the hard state activity of the source was detected with the IBIS/IS GRI imager and the JEM- X 
telescope on board the INTEGRAL satellite on 7—8 February, 2011 (jCapitanio et al.ll201ll ). They 
described the combined ISGRI and JEMX2 spectrum (5.0-300.0 keV) with a cut-off power law 



model with the high energy cutoff of llOl^^ keV. 

From the beginning of the RXTE observations on 9 February 2011, several rich features 
and characteristics in timing and spectral domain of the source have been noticed. QPOs at 
0.094 ± 0.003 Hz with a quality factor of ~ 9.1, and at 0.105+o;oo3 Hz with a quality factor 
of ^ 7.4 etc. were dete cted in the RXTE/PCA observations on 9, 10, 11 and 12 February, 2011 



47 ( Rodriguez et al.ll2011all. Apart from this, fast QPO evolution (|Shaposhnikovll201ll ). 1 mHz QPOs 



Altamirano et al.ll2011al ) and evolving power density spectr um with multiple peaks (IPahari et al 



201 Ibl ) were also noticed. High frequency QPOs at 66 Hz (lAltamirano &: Bellonil |2012|) were de 



50 tected with the significance of 8.5c7. Using Chandra spectroscopy of IGR J17091-3624, iKing et al 



(|2012l ) showed that accretion disk wind in soft, disk-dominated states, in apparent anti-correlation 
with relativistic jets in low hard state. This disk wind is observed to be highly ionized, dense, and 
has typical blueshift of ~9300 km/s or less projected along our line of sight. Lately, a quasi-regular 
oscillation in the lightcurve of IGR J 1709 1-3624 was noticed which is very simil ar to the '/o' class 



55 variab ility seen only in the lightcurve o f galactic microquasar GRS 1915+105 (jAltamirano et al 



2011bl ). Later, lAltamirano et al.l (j201ld ) claimed few more variabilities similar to those found in 
GRS 1915+105. Two new variabilities which are unique to IGR J17091-3624 are reported by 



58 (jPahari et al.l l2011al ) . 



GRS 1915+105 is a Galactic microquasar which has shown spectacular X-ray variability since 
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its discovery in 199 2 (ICastro-Tirado et al. I I1994I ). Out of 12 classes observed in the X-ray variabilities 
of GRS 1915+105 (jBelloni et al.ll200d ). the class shows regular, repetitive bursts in the lightcurve 
with recurrence time between 40s and 110s. These recursive burst structures have been observed 



63 for extended duration a nd have been discussed by several authors (jTaam et alj 119971 : iPaul et al 



19981 : lYadav et al.1 1 19991 ) . Using the data from the Indian X - ray A stronomy Experiment (IXAE) 



onboard the Indian Satellite IRS-P3 and RXTE, Yadav et al. (j 19991 ) found that these bursts which 
can be categorized as irregular, quasi-regular a nd regular burs t s have a slow exponential rise and 
sharp linear decay. Using the BeppoSAX data, iMassaro et al.l (|20ld ) concluded that the process 
responsible for the pulses produced strongest emission between 3.0 and 10.0 keV, and the emission 
at the rising phase of the pulse dominated the lower and higher energy. However, this hypothesis is 
not tested. Hence origin of 'p' class activity in GRS 1915+105 is poorly understood. The similarity 
of the p X-ray class of GRS 1915+105 with the variability observed in IGR J17091-3624 motivated 
us to perform a detailed analysis of the current outburst of IGR J17091-3624. 

In the present work, we have considered all observations of this source between 09 February 
2011 and 10 April 2011 with the RXTE/FCA and the SWIFT/XRT. We analysed both timing 
and spectral properties of the source after confirming that the contamination from the nearby 
source IGR J17098-3628 was usually very small. In section 2 we describe our analysis methods. 
Spectral analysis of the SWIFT/XI{T show the steepening of powerlaw index when the X-ray flux 
(0.3-10.0 keV) increased by a factor of ~ 4.5. A gradual and prominent change in the variability 
was detected which led to the class transition of the source. The nature of the compact object 
(in X-ray binaries) plays an important role in selecting spectral model parameters. However, na- 
ture of compact object in IGR J17091-3624 is not well determined. Spectra of black hole X-ray 
binaries (BHXBs) are generally described by a soft X-ray emission from the disk and hard power- 
law emission from cor ona/base of jet. Using combined RXTE/FCA and INTEGRAL/IBIS data. 



Capitanio et al.l (|2006l ) found that a single thermal comptonization component (compTT in XSPEC) 



or saturated comptonization (compST in Xspec) or the combination of soft multicolor disk black- 
body and hard powerlaw (diskbb+powerlaw) are good description of the X-ray spectra of IGR 
J17091-3624. However, they found the electron temperature in thermal comptonization model 
was ~ 20 k eV, which was less than typ ical electron temperature of black hole X-ray binary (~ 
50-100 keV) (jZdziarski fc Gierlinskil |2004| ) . The hard X-ray outburst in this system and the large 
flux of powerlaw spectral component over soft thermal component are commonly seen in BHXB 
systems. The preliminary spectral data during the recent outburst in 2011 showed that a soft spec- 
tral component is nec essary to describe the spectrum along with the hard powerlaw component 



93 (Del Santo et al 



2011 



Section 3 summarizes all results regarding hard state and the hard inter- 
mediate state, fast source evolution and transition to a 'p'-like class via intermediate variability 
state, observed in IGR J17091-3624. Finally we have compared the timing and spectral properties 
of the '/o'-like class with the 'p' class characteristics seen in GRS 1915+105 in section 4. We give 
our conclusion in section 5. 
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98 2. Observations and Data Analysis 

99 We analyzed all 34 pointed observations of IGR J17091-3624 made with RXTE/PCA between 

100 9 February 2011 (MJD 55601) and 10 April 2011 (MJD 55661). We used archival data for 28 

101 observations from 9 February 2011 to 29 March 2011, and real time data for last six observations. 

102 Due to insufficient amount of data and inferior quality, we were unable t o extract the backgroun d 



103 data file for all of the High-Energy X-ray Timing Experiment (HEXTE: iRothschild et al.l (|l998l )) 

104 observations. Hence, we used only the RXTE/PCA standard-2 science array data and event mode 

105 good xenon data for our analysis. Observation details are given in Table 1. We also analyzed 

106 three RXTE/FCA observations of GRS 1915-^105 on 26 May 1997, 03 June 1997 and 22 June 

107 1997 when the sources was in 'p' class. RXTE/PCA consists of an array of five proportional 

108 counters (PCUO-4) filled w ith xenon gas, with total effective area of ~ 6500 cm^, operating in 



109 2.0-60.0 keV energy range (jjahoda et al.lll996l ). We used data observed with all xenon layers of 



no PCU2. During the start or the stop of any PCU unit, many observations showed large count rate 

111 fluctuations due to the instrument in the lightcurve. To take care of that, we first created xte filter 

112 file for each PGA observation and calculated good time interval (GTl) from the exposure time 

113 by applying all standard filtering criteria including breakdown events. Galculated values of good 

114 time interval are provided in Table 1. Using GTI values of individual observation, we extracted 

115 background subtracted light-curve from good xenon data having resolution of 125 micro-sec. As the 

116 variability time scale found in light curve is larger than the actual resolution, we used Is re-binning. 

117 Observations with the background subtracted PGA intensity > 15 cts s~^PGU~^ were considered. 

iiB We also analyzed 27 archival data sets of SWIFT/X.P!T between 03 February 2011 and 30 March 

119 2011. The XRT instrument on board S'VF/FT satellite has the effective area is 110 cm ^ at 1.5 keV, 



120 the p osition accuracy is ~ 5 arcsec and is operating in 0.2-10.0 keV energy range (jBurrows et al 



121 I2OO5I ). The typical spectral resolution of this instrument is ~ 140 eV at 6.0 keV. To avoid pile- 

122 up problem while observing luminous X-ray sources, all observations (except the first one) were 

123 performed with Windowed Tir ning (WT) rnode w hich has time resolution of 1.8 ms and flux limit 



124 of ~ 0.6 mGrab (~ 200 cts/s) (jMineo et al.ll2006l ). For our analysis, we used level 2 cleaned event 



125 data which had been extracted from level 1 calibrated data after applying several screening criteria 

126 on specific parameters like GGD Temperature, Sun Angle, Elevation Angle etc. For all WT mode 

127 data, we assigned grade 0-2 and selected good event with STATUS==bO. For each observation, we 

128 first extracted the image. Then following the standard procedure for selecting source region as well 

129 background region, we extracted source and background lightcurves and spectra separately using 

130 xselect v2 . 4b in FTOOLS 6.10. 

131 Top panel of Fig. 1 shows that the background-subtracted average RXTE/PCA count rate 

132 between 2.0-60.0 keV. We found sudden drop in RXTE/PCA count rate from ~ 600 cts/s/PGU to ~ 

133 190 cts/s/PGU after first 14 days of observations. Surprisingly, at the same time, the background- 

134 subtracted average SWIFT/XKT count rate between 0.3-10.0 keV showed an increase by a factor 

135 of ~ 4. Detailed analysis showed that during first 14 days of observation, IGR J17091-3624 was 

136 within the field of view of persistently luminous Z-type neutron star low mass X-ray binary GX 
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137 349+2 (17''05'"44.49^ (R.A. (J2000.0)), -36°25.38' (decimation (J2000.0))). As a consequence, very 

138 high PCA count rate was observed in first 14 days of observation compared to the rest. Therefore, 

139 in our present work, we considered PCA observations from 23 February 2011 onwards (observations 

140 falls right to the straight line in the top panel of Fig. 1) to avoid the contamination effect. Further 

141 to check possible contamination from IGR J17098-3628, a transient black hole X-ray binary located 

142 9.6 away from IGR J17091-3624, we analyzed background-subtracted image from SWIFT/XBT in 

143 the energy range 0.3-10.0 keV (Fig. 2). The position of IGR J17098-3628 (showed by a circle) is very 

144 faint compared to the position of IGR 224 J17091-3624 (showed by a square). For a quantitative 

145 measurement, we calculate the count rate separately from the position of IGR J17091-3624 and 

146 from the position of IGR J17098-3628. For the energy range 0.3-10.0 keV, their values were 2.023 

147 cts/s and 0.213 cts/s on 3 February 2011. Hence, IGR J17098-3628 has the intensity ~ 10% of 
14B IGR J17091-3624. Since, RXTE/FCA could not resolve these sources, it is expected that ~ 10% 

149 contamination effect would be present in RXTE/FCA data. 

150 2.1. Timing Analysis 

151 For each PCA observation we created a deadtime corrected and white noise subtracted power 

152 density spectrum (PDS) to study different variability features and to track the nature of change 

153 occurred in the source with time. The dead time corrections are first applied to the lightcurves. 

154 For GoodXenon data, the dead time per event is approximately 10 /is. For details of the procedure, 

155 please check RXTE cookbool<0. These corrected lightcurves were used to create PDS in 2.0-60.0 

156 keV energy range. To improve the signal-to- noise ratio in the PDS, we used the geometric re- 

157 binning of frequency bins of the factor of 1.0^. The expected white noise was subtracted from the 

158 PDS^ and they were normalized such that their square root of integral over the range of frequencies 

159 provide fractional rms variability. For same values of bin time and geometric re-binning we also 

160 produced PDS from the SWIFT/^KT data using background subtracted light-curve between 0.2- 

161 10.0 keV energy range. We used same method of normalization. Since one could not compare 

162 the 2.0-60.0 keV PDS from the PCA data with 0.2-10.0 keV PDS from SWIFT/XRT data, hence 

163 we extracted SWIFT/XI{T PDS for observation dates for which only SWIFT/XI{T data were 

164 available. Several QPO-like features were noticed in the PDS of different RXTE/FCA observations. 

165 Their approximate characteristics are given in Table 1. Most of the S'VFIFT/XRT PDS were noise 

166 dominated or without any feature. 



^http: / /heasarc. nasa.gov / docs /xte / recipes / pca_deadtime.html 
^http: / /heasarc. nasa.gov / xanadu / xronos /help /powspec.html 
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167 2.2. Spectral Evolution 

16B The color-color diagram (CD), the hard color vs. soft color plot and the hardness intensity 

169 diagram (HID) were created using the PCA Good Xenon data with 1.0 s bin time. The soft color 

170 was computed by dividing the background subtracted PCA count rates in the energy band 5.0-12.0 

171 keV by that in the energy band 2.0-5.0 keV. Similarly, the hard color was defined as the ratio of 

172 the background subtracted count rate in the energy band 12.0-60.0 keV and 2.0-5.0 keV. X-ray 

173 intensity is the background subtracted PCA count rate in 2.0—60.0 keV energy range. Since the 

174 channel gain varies from epoch to epoch, we calculated absolute channel ranges corresponding to 

175 each of the energy ranges used to calculate soft color, hard color and X-ray intensity for the present 

176 epoch. For each RXTE/PCA observation, mean count rate (top panel of Fig. 1), and mean value of 

177 hard color and soft color have been calculated for a given set of data and their standard deviation 

178 values were provided as error (provided in table 1). 

179 2.3. Spectral Fitting 

180 In case of RXTE/PCA, we performed the time segmented spectral analysis using GoodXenon 

181 data. The details of the analysis is given in section 3.4.1. In case of SWIFT/XI\T data, we first 

182 extracted raw ID image of the source. Then choosing suitable source region and background region, 

183 we extracted source spectrum and background spectrum separately from the cleaned and pointed 

184 WT mode event file. We used the latest standard response matrix file (rmf) and auxiliary response 

185 file (arf ) for WT mode data between 0-2 grade values. We used XSPEC v 12 . 6 . for spectral fitting 

186 of RXTE/PCA and SWIFT/'XJ{T data. Sudden drop in spectral counts with poor signal-to-noise 

187 ratio have been observed for lower energy channels of PCA. Similarly very high energy channels 

188 are also affected by high background count rate. Hence to maintain uniformity, we used ^ 2.5-28.0 

189 keV energy spectrum for our spectral fitting. 

190 We have tried different single component models like diskbb or compTT or powerlaw as well 

191 as dual component models like diskbb+powerlaw, diskbb+compTT or compTT+powerlaw to fit the 

192 source spectra. The reason for selecting and testing different models are greatly discussed in [H We 

193 found that the single SWIFT/XBT spectra could be fitted well with a combination of diskbb and/or 

194 powerlaw model. Fitting with the rest of the model combination yielded unacceptably large reduced 

195 x^- We used a small Gaussian component at 6.4 keV in order to improve the reduced iii some 

196 observations. To account the effect of absorption by the neutral Hydrogen, all model components 

197 were multiplied by a photo-electric absorption model. The Galactic absorption column density at 

198 the direction of 17''11'"12.5'' (R.A. (J2000.0)), -36°24.6' (declination (J2000.0)) was calculated as 

199 0.79 X 10^^ cm^^. Again from spectral analysis, we found that for all observations, absorption 

200 column density varies between 0.82 it 0.2 and 0.99 it 0.3. Thus we kept this value fixed at 0.9 x 

201 lO^'^ cm~^ for all spectral fitting. The unfolded fitted spectra between 2.5-28.0 keV are shown in 

202 bottom left panel of Fig. 3 along with its model component and residual. Table 2 shows the fitted 
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203 values of different model component parameters. 



3. Results 

After determining the part of observations which is useful for our present work, we performed 
the detailed analysis of the combined RXTE/FCA and SWIFT/XRT data. SWIFT/XRT observa- 
tions of IGR J17091-3624 showed a transition of the source from the low count rate hard state to 
the high count rate hard state (see top panel of Fig. 1) and a gradual transition in the variability, 
i.e., from an irregular non variability state, to a regular, repetitive and large variability state. This 
new variab ility looked very similar to the unique variability observed during the 'p' class of GRS 



211 1915+105 (jAltamirano et al.ll2011bl ). This motivated us to perform a comparative study between 

212 these sources. For better distinction of above states, we performed the count rate fractional rms 

213 analysis of different observations (by dividing the square root of the variance with the mean count 

214 rate) from RXTE/PCA which is shown in the bottom panel of Fig. 1. Observations which have 

215 fractional rms values less than 3.2 are found to be in hard intermediate state and fractional rms 

216 more than 3.2 are found to be in variable state. A few observations showed the transition between 

217 hard intermediate state and variable state, and hence described as 'Intermediate variable state 

218 (IVS)'. The details have been discussed below. 



219 3.1. Hard State and Hard intermediate state 

220 From 14 February 2011 to 24 February 2011, a gradual transition in the SWIFT/XRT flux 

221 between 0.3-10.0 keV was observed in IGR J17091-3624. Initially till 14 February 2011, average 

222 SWIFT/XRT count rate was ~ 10 cts/s. Then it gradually increased to ~ 45 cts/s on 24 February 

223 2011. To check whether this rise in flux led to any state transition, we analyzed the SWIFT/XRT 

224 data of 10 February 2011 when the flux is low and 04 March 2011 when the flux is high. Background- 

225 subtracted light curve for both observations did not show any specific change in variability (Top left 

226 panel of Fig. 3). Secondly, the Leahy normalized power density spectra of both observations show 

227 similar spectral shape and the total integrated power between 0.1 Hz and 10.0 Hz are also close to 

228 each other (top right panel of Fig. 3). Hence timing properties of the source remain unaffected by 

229 the increase of the count rate. 

230 Now to check the spectral behavior of the source during low-flux to the high-flux transition, 

231 we did energy spectral analysis of both observations, i.e., 10 February 2011 and 04 March 2011 

232 respectively. We tried different model components the description of which is given thoroughly in 

233 § 12.31 While fitting with 1% systematic error, different models described in § 12.31 failed to describe 

234 the spectra except diskbb+powerlaw model. All of them (like diskbb, compTT, diskbb+compTT 

235 etc.)fitted the spectra with unreasonably high reduced value (>1.5) while diskbb+powerlaw 

236 model gives acceptable fit with reduced ~0.85. This is not so surprizing, since most of the BHXB 
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spect ra are well described by the diskbb+powerlaw model combinations (jRemillard &: McClintock 
20061 ). However, we found that disk component is not significant while fitting the spectrum observed 
on 10 February 2011. A single powerlaw can adequately describe all low flux spectra between 03 
February 2011 and 22 February 2011. Bottom panels of Fig. 3 show the unfolded energy spectra, 
fitted model components and residuals of observations during low flux on 10 February 2011 (left) 
and during high flux on 04 March 2011 (right). Table 2 show the parameter values of fitted model 
components during the low flux state and the high flux state respectively. From Table 2, it is clear 
that all low flux state observations have photon index value less than 1.8 and no disk component 
is present at the detectable limit o i SWI FT/XI{T. These meet the criteria of a Hard state as 
defined bv iRemillard McClintockl (j2006l ). On the other hand, all high flux state observations 
have powerlaw indices between 2.1 and 2.4. Now spectra with powerlaw index between 2.1 (upper 
limit for hard state) and 2.4 (l ower limit of the Steep Power La. w (SPL) state) are defined as a Hard 
intermediate state (HIMS) bv IRemillard &: McClintockl (j2006l ). The increase in powerlaw photon 
index is first observed on 22 February 2011 using SWIFT/XWT (See Table 2). Later we confirm this 
by performing RXTE/PCA spectral analysis of all observations under consideration. We found that 
all observations with low flux belong to the hard state while all observatio ns with high flux belong 
to the HIMS. Our results agree well with spectral parameters obtained by lCapitanio et al.l (120121 ) 
using the joint spectral analysis of INTEGRAL/lBlS and SWIFT/XRT data. The relatively smah 
difference between our results may be due to their (1) use cutoff powerlaw (where cutoff energy >70 
keV, see Table 1 of ICapitanio et al.l (|2012l )) (2) different absorption column density (1.1 it 0.3). 



Using RXTE/FCA rms normalized power density spectra between 2 mHz and 12 Hz, we found 
a significant evolution in the QPO frequencies in different observations. Table 1 show that from 
26 February 2011 to 6 March 2011, QPO frequency increased from ~ 3.81 Hz to ~ 5.53 Hz. Then 
it decreased to ~ 5.17 Hz on 10 March 2011. No high frequency QPOs were significantly detected 
upto 100 Hz. 



262 3.2. Intermediate variable state 

263 During the continuous monitoring of the PCA light curve of IGR J17091-3624, we found a 

264 change in the variability pattern from 02 March 2011 to 19 March 2011 (see the top left panel of 

265 Fig. 4) in the energy range 2.0-60.0 keV. On 02 March 2011, the source was in hard intermediate 

266 state when the light curve was less variable and very similar to the typical light curves seen in 

267 canonical BHXBs. On 19 March 2011, the light curve showed a very regular repetitive pattern and 

268 high variability (peak count rate ~ 3—4 times higher than the count rate at the persistent level). 

269 We denoted it by variable state. Observations and average count rate that belonged to the HIMS 

270 and variable state were separated by a dashed line in Fig 1. 

271 This transition occurred through some semi-oscillatory intermediate stage on 14 March 2011. 

272 The random noisy pattern in the lightcurve of 02 March 2011 seems to take an oscillatory shape on 

273 1 4 March 2011. Although these intermediate phases are transient but it has been observed couple of 
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274 times like on 23 & 24 March 2011, when the source oscillates between the HIMS and variable state. 

275 The rms normalized PDS during the intermediate variable state (IVS) showed QPOs in 10—20 mHz 

276 frequency range which increased to 20—50 mHz range during variable state. HIMS PDS do not 

277 show any such mHz QPO which remarkably distinguish the variable state from the HIMS (Table 1). 

278 Interestingly, the HIMS PDS and the IVS PDS share equal source power continuum except a strong 

279 QPO at 4.18 Hz during the HIMS whereas the variable state PDS shows very high power with mHz 

280 QPO below 0.5 Hz and it drops to noise power continuum before 1 Hz. Fractional rms amplitude 

281 of the PDS was found to be 11.6% which is intermediate between 5.8% (HIMS) to 21.4% (variable 

282 state) (top right panel of Fig. 4). To trace any change in spectral state during this transition, 

283 we plotted in Fig. 4 the HID (bottom left panel) and the CD (bottom right panel). In the HID, 

284 the intermediate variable state shows slightly higher flux and larger hardness variability. Prom the 

285 HID, it is seen that the lower limit of the hard color value in the intermediate phase is less than the 

286 lower limit of the HIMS. In the CD, the distribution pattern slightly moved towards bottom left 

287 corner. However, in both cases, intermediate variable state lie between HIMS and variable state. 

288 These results show significant change in the timing properties of the source but not much change in 

289 spectral properties. To check this, we performed energy spectral analysis of Intermediate variable 

290 state (IVS) on 14 March 2011 between 2.5—30.0 keV. The spectrum can be well fitted by diskbb 

291 and powerlaw model which showed parameter values similar to the HIMS. 

292 3.3. Variable state 

293 The new variable state remained stable for almost 26 days. The recurrent time scale of this 

294 variability was found to change with time. Fig. 5 showed the evolution of the source via five 

295 RXTE/PCA hghtcurves from 19 March 2011 to 10 April 2011. On 19 March 2011, the variability 

296 timescale of recursive bursts was ~ 41 sec and the burst peak count rate was ~ 2.92 times more than 

297 the count rate at minimum level. These bursts have simple profile with slow-rise and fast-decay. 

298 As these bursts evolve, flaring frequency increases. On 10 April 2011, the recurrence time scale of 

299 the variability became ~ 21 sec. and the burst became more structured. The ratio of burst peak 

300 count rate to the persistent level count rate became ^ 3.38. Careful observation of burst profiles 

301 revealed two unique burst structures. (1) Burst profile seen on 31 March, 2011. They have a single 

302 peak and the count rate rises slowly (~ 30 sec) to the peak of the burst. (2) Burst profile seen on 

303 03 April, 2011, these bursts are characterized by double peak and show faster rise in count rate (~ 

304 15-20 sec). Their peak-to-persistent level count rate ratio is ~2.84. This behavior and variation in 

305 burst profile structures of the source has so far been seen only in the 'p' class of GRS 1915-1-105. 

306 Prom the bottom panel of Fig. 1 and Fig. 5, it is clear that variabilities becomes stronger (i.e., 

307 higher amplitude of oscillation) and faster (i.e., smaller recurrence time-scale of the burst) with 

308 time. From the PDS, we found dual QPO-like features in few observations. For example, QPO 

309 features were found at ~ 26.6 mHz and ~ 51.8 mHz on 22 March 2011, and at ~ 24.8 mHz and 

310 ~49.6 mHz on 29 March 2011, at ~ 37.2 mHz and ~ 75.3 mHz were found on 03 April 2011, at ~ 

311 46.5 mHz and ~ 94.0 mHz were found on 10 April 2011 etc.) (see Table 1). Multiple QPOs were 
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312 observed at ~ 24.3 mHz, 45.2 mHz & 70.4 mHz on 30 March 2011 and at ~ 28.1 mHz, 54.5 mHz 

313 and 84.6 mHz on 31 March 2011 (see Table 1). Noise components with variable strength were 

314 detected in different observations. We also detected QPOs around 8—10 Hz during this state (see 

315 Table 1 and upper middle left panel of Figure 9). 

316 In the following section, we explored other properties of the variable state, and compare them 

317 with those of GRS 1915+105. 



31B 3.3.1. Properties of variable state in comparison to GRS 1915+105 

319 As the variable state in IGR J17091-3624 looks similar to the 'p' class of GRS 1915+105, we 

320 make a comparative study between them. 

321 All panels of Fig. 6 showed the 2.0-60.0 keV light curve in both sources. To track the change 

322 in variability pattern, we took two observations of variable state in IGR J17091-3624 which show 

323 transition from the slow & low flux variability (31 March 11) (top left panel of Fig. 6) to the 

324 fast &: high flux variability (10 April, 11) (bottom left panel of Fig. 6). Similar transition from 

325 the slow & low flux variability (top right panel of Fig. 6) to the fast &: high flux variability 

326 (bottom right panel of Fig. 6) was observed in GRS 1915+105 on 26 May 1997 and 22 June 1997 

327 respectively. However, it is interesting to observe that that the peak count of the bursts in GRS 
32B 1915+105 increased significantly while transiting from the slow variability to the fast variability 

329 (~ 4000 cts/s/PCU to ~ 5000 cts/s/PCU) while it increased slightly in IGR J17091-3624 during 

330 the variable state. This variability in IGR J17091-3624, although repetitive, was irregular, less 

331 structured and more quasi-periodic compared to the 'p' class of GRS 1915+105. 

332 Previously a strong anti-correlation between the variability pattern in lightcurve and hardness 

333 ratio (defined as the ratio of count rate betwe en 12.0-60.0 keV and 2.0-12.0 keV) was found in the 



334 class of GRS 1915+105 (jYadav et al.lll999l ). This anti-correlation was also found in the variable 

335 state of IGR J17091-3624. Fig. 7 shows the comparative result of anti-correlation found in IGR 

336 J17091-3624 on 31 March 11 (top and bottom left panels) and in GRS 1915+105 on 22 June 97 

337 (top and bottom right panels). It is clear that the hardness ratio is higher in IGR J17091-3624 
33B than that observed in GRS 1915+105. 

339 To study the nature of the bursts in details, we studied the rise and the decay profiles for both 

340 sources. Left panel of Fig. 8 showed the combined rise profiles of several bursts and combined decay 

341 profiles of same bursts in IGR J 1709 1-3624 on 31 March 11. To study the average behavior of rise 

342 profiles, we normalized the starting time of each profile to sec. Similarly, to study the average 

343 behavior of decay profiles, we normalized the starting time of each profile to 15 sec. We repeated 

344 the same procedure for the observation on 22 June 97 of GRS 1915+105 (shown in the right panel of 

345 Fig. 8) except for the fact that the staring time of each decay profile was normalized to 70 sec. We 

346 fit the combined rise profile with exponential function frise(t) = ae*/*''''"= (red continuous line), the 

347 combined decay profile with exponential decay function idecayi^) = ae~*/*-f°" (blue continuous line) 
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348 and all profiles with straight line (red and blue dotted line). In case of GRS 1915+105, exponential 

349 function fits better than straight line in both rise and decay profile with the significance of 1—1.94 

350 X 10"^'' and 1—1.31 x 10"'^^ respectively. In IGR J17091-3624, the corresponding significance 

351 values were 0.312 and 1—1.99 x 10~^ respectively. The slope of the exponential function, for both 

352 rise and decay, were found similar in both sources (see Table 3). This indicate that the origin of 

353 burst structure might be similar in both sources. 

354 Top panels of Fig. 9 showed that the low frequency (0.1—10.0 Hz) power, QPO strength and 

355 the nature of noise components are similar for both sources. White noise subtraction was not 

356 applied only for PDS shown in the top panels of Fig. 9. High frequency power disappeared above 

357 1 Hz in case of IGR J17091-3624(top left panel of Fig. 9) but the powerlaw like continuum in 
35B GRS 1915+105 continued at higher frequencies along with broad QPO like feature at 8—10 Hz 

359 (top right panel of Fig. 9). Harmonics at mHz frequency are seen in both sources. In case of IGR 

360 J17091-3624, we have found 8-10 Hz QPO during 'p' like X-ray activity on 25-03-11, 27-03-11, 

361 31-03-11 and 10-04-11 (see Table 1). Upper middle panel of Fig. 9 show 8-10 Hz QPOs for both 

362 sources. It is notable that the strength of 8—10 Hz QPOs in IGR J17091-3624 are lower than that 

363 observed in GRS 1915+105. Lower middle panels of Fig. 9 show the comparative study of the 

364 HIDs between two sources which shows similar general trend but slightly different range in hard 

365 color for both sources. For the quantitative measurement of the variability evolution, we considered 

366 the PDS of few observations where the first and the last observation were 26 days apart for both 

367 sources. During all observations under consideration both sources were in 'p'-like class. Bottom left 

368 panel of Fig. 9 showed that the fundamental peak in the PDS of IGR J 1709 1-3624 shifted from ~ 

369 12.3 mHz (on 14 March 2011) to ~ 46.8 mHz (on 10 April 2011) in 26 days. While in case of GRS 

370 1915+105, the shift of the fundamental peak in the PDS is from ~ 9.7 mHz (on 26 May 97) to ~ 

371 18.5 mHz (on 22 June 97). Thus the evolution was much faster in IGR J17091-3624 (the change in 

372 QPO frequency is ~ 34.5 mHz in 26 days) than GRS 1915+105 (the change in frequency is ~ 8.8 

373 mHz in 26 days). Results on comparative study of both sources are summarized in Table 3. 



374 
375 
376 
377 
378 
379 
380 
381 
382 
383 
384 
385 
386 



During 2011 outburst, IGR J17091-3624 evolved from the hard state (upto 22 February, 2011) 
to the HIMS state and finally to the variable s tate ( similar to p class in GRS 1915+105) as X-ray 
flux increases (see Table 1 of ICapitanio et al.l (120121 ) & our Table 2). Detection of accretion disk 
wind as well as quenching of radio emis sion during disk dominated state suggests that IGR J17091- 
3624 is approaching the high soft state (jKing et al.ll2012l : [Rodriguez et al.ll2011bl ). The HIMS state 
is characterized by by 3.80—7.08 Hz QPO as well as powerlaw index in the range 2.1—2.4 (Table 1 
& 2). During 1997, GRS 1915+105 showed p class activity during May-June, 1997 when it evolved 
from the HIMS state (bright low hard state observed on 26 March, 19 97 &: 8 May, 1997) to the high 



soft s tate during July- August, 1997 (20 July, 1997 & 19 August, 1997) (jRao et al.ll200d : lYadav et al 



19991 ). In GR S 1915+105, HIMS shows QPOs in the range 3.5-4.2 Hz and powerlaw index in the 



range 2.1— 2.5 ()Rao et al 



200C 



Muno et al.lll999l ). GRS 1915+105 is never seen in the typical low 



hard state (F < 1.6) since its discovery in 1992 but frequently seen in very high state/luminous SPL 
state where transient radio jets are produced (|Mirabel et al.lll998l ; lEikenberry et al.lll998l ; lYadav 
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20061 ). It may be noted that IGR J17091-3624 is seen in the hard state but not seen in very high 



state/luminous SPL state. Both these sources show p class activities when evolve from the HIMS 
state. Other BHXBs which show transient radi o jets evolve from the quiescen t state/low hard 
state to the very high state/luminous SPL state (jRemillard &: McClintockl 120061 ). Therefore, IGR 
J17091-3624 is providing missing link between GRS 1915+105 and other BHXBs. 



392 4. Discussion 

393 We have studied the evolution of the recent X-ray activity in IGR J17091-3624 from 3 February 

394 2011 to 10 April 2011. Initially a hard to hard intermediate state (HIMS) transition was detected 

395 by the SWIFT /^KKT within 0.3—10.0 keV energy range. The hard state is mainly characterized 

396 by the low count rate, power-law dominated energy spectrum with the powerlaw index less than 

397 2.1 and the power density spectrum with very high rms. The HIMS is characterized by the high, 

398 nonvariable count rate, powerlaw index between 2.1 and 2.4 and both disk component and powerlaw 

399 component have significant contribution to the total flux. Few days later, a transition from the 

400 HIMS to a regular, repetitive and highly variable state occurred. A few observations between 

401 variable state and the HIMS showed low intensity quasi periodic variabilities which we have termed 

402 as intermediate variable state. Spectroscopically they are similar to the HIMS. After that, nearly 26 

403 days, the source showed regular, structured, repetitive variabilities similar to the p-class, observed 

404 only from GRS 1915+105, a Galactic microquasar harboring a black hole of mass 14.0+4 M©. 

405 In addition to the structural similarity with GRS 1915+105, we found that both sources during 

406 p class showed similar evolution in the burst structure. With time, burst frequency increases in 

407 both sources as well as the ratio of the peak flux to the minimum level flux also increases. In both 

408 cases, average rise profile of the burst shows slow exponential rise while average decay profile shows 

409 fast decay. Interestingly, we found strong anti-correlation between hardness value and the X-ray 

410 flux during variabilities in both sources. The evolution of the burst structure in the HID of both 

411 sources is identical. PDS below 1 Hz, showed similar structure. Both sources have strong mHz 

412 QPOs sitting on the top of the white noise and both of them showed steep red noise in the PDS 

413 along with 8-10 Hz QPOs. 

414 In other BHXBs, transition from the canonical low hard state to the steep power law state 

415 via hard intermediate state is commonly observed phenomenon. IGR J17091-3624 also showed 

416 transition from the hard state to the HIMS. Spectral properties of the hard state and the HIMS 

417 of IGR J17091-3624 also matches with that of other BHXBs. However, GRS 1915+105 is usually 

418 seen in the HIMS state (but never seen in the typical low hard state) and it moves to the very high 

419 state/luminous SPL state via high soft state. Again, IGR J17091-3624 showed regular repetitive 

420 variabilities similar to the p class variability observed in GRS 1915+105. This kind of variability has 

421 never been observed in any other BHXBs. Thus IGR J17091-3624 showed few properties similar to 

422 GRS 1915+105 as well as similar to other BHXBs. Hence it is behaving like a missing link between 

423 GRS 1915+105 and other BHXBs. 
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424 Although similar observations do not necessarily imply that both sources are of similar type 

425 but we found some evidences which tentatively indicate that the central object in IGR J17091-3624 

426 system may be a black hole candidate. One such indication is the detection of type-B QPO at 

427 least four times and the detection of type-C QPOs several times. Transient-type Type-B QPOs 

428 were seen in few BHXB systems in the frequency range 5—6 Hz around state transition and it has 

429 low fractional rms a mplitude (2— 4 %) a nd low quality factor ~ 6 and associated with weak red 

430 noise (for details, see ISoleri et al.l ([2003)). We detected four QPOs at 5.38, 5.53, 5.39 and 5.17 

431 Hz during the HIMS. Since no PC A data are available during the hard state, hence QPO analysis 

432 during hard to HIMS transition is not possible. They have average fractional rms amplitude about 

433 2—3% and quality factor around 4—5 (for details see Table 1). All of them are associated with 

434 weak red noise. Thus they satisfied all conditions of type-B QPO. Type-C QPOs are commonly 

435 found in BHXB systems and typically found in the frequency range 0.1—15 Hz. They have high 

436 rms amplit ude (7—12 % ) , high quality factor (~ 10) and associated with strong flat top noise (for 

437 details, see ISoleri et al.l (|2008l )). We found several of them in our analysis and type-A QPO was 

438 found once at ~ 7.06 Hz (see Table 1). Non-detection of upper and lower kHz QPOs in this source 

439 indirectly decreases the ch ance for the source of b eing a neutron star. Again the break frequency 

440 in the PDS of this source (jLutovinov et al.ll2003l ) further decreases the chance. No type-I X-ray 

441 bursts were found in any observations within two months. These bursts unique characteristics of 

442 neutron star LMXBs. This gave an indirect evidence that the source might not be a neutron star. 

443 However, we found few dissimilarities between the variable state of IGR J17091-3624 and the 

444 '/)' class of GRS 1915-1-105. They are listed below: 



445 
446 
447 
448 



The c olor-color diagram (CD) of the 'p' class showed a loop-like pattern (see iBelloni et al 
([20001)) while the CD of the variable state in IGR J17091-3624 showed a patchy pattern (see 
middle right panel of Fig. 4). The X-ray flux in IGR J17091-3624 is lower than that observed 
in GRS 1915-^105. 



449 

450 
451 
452 
453 
454 
455 



An interesting difference between two systems is that the hardness ratio during p class activity 
is higher (~ 2 times) in IGR J17091-3624 than that observed in GRS 1915-M05. Adjustment 
of black hole mass and its distance may not be sufficient to explain this difference. Peak-to- 
minimum level count rate ratio for IGR J17091-3624 is ~ 2.9—3.8. while the value of this ratio 
is ~ 3.4—4.3 in case of GRS 1915-1-105. The absorption column density and contamination 
from close by source are other parameters which may affect hardness ratio. However, the 
chance of contamination from close by source is low. 



456 
457 
458 



Typical low hard state observations where disk emission is not visible, are commonly seen in 
canonical BHXBs but have never been seen in GRS 1915-1-105. However, IGR J17091-3624 
showed low hard state from 9 March 2011 to 22 March 2011 (see Table 2). 



459 
460 



4. 



A systematic study of t he occurrence of different classes during the outburst of GRS 1915-1-105 
([Yadav &: Paharill201ll ) showed that the p class variability eventually arise from the q or a; 
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464 



479 
480 



class (jPahari &: Palll20ld ). These are different X-ray classes with different timescale and flux. 



462 However, in IGR J17091-3624, the p class like activity arises from the HIMS state. 

463 5. The variable state was found stable quasi-continuously for ~ 26 days in IGR J17091-3624, 
while it was found to be stable for 10—30 days in the GRS 1915+105 alt hough it was not 



contin uous. These observations, in both sources, proved the conclusion of lYadav Pahari 



(j201ll ) that class variability represented very stable X-ray class. The possible reason for 
prolonged stability of IGR J17091-3624 in the variable state compared to GRS 1915+105 was 
that to trigger the transition from one variability to another, the disk should cross a critical 
threshold mass accretion rate. The timescale for attaining this critical mass accretion rate 
might be longer which prolonged the occurrence of 'p' class in IGR J17091-3624. 



Using the radio data, iRodriguez et al.l (1201 Ibl ) found the distance of IGR J17091-3624 to be 



472 10—17 kpc assuming black hole mass to be 10 M©. We used the distance of IGR J17091-3624 to be 

473 14+3 kpc and the compact object mass to be 10 M0 to calculate the mass accretion rate of IGR 

474 J17091-3624 relative to the GRS 1915+105. In case of GRS 1915+105, we assume the mass, the 

475 distance and the disk inclination angle with respect to observer's line of sight to be 14.0+4 Mq, 

476 12+2 kpc and 70° respectively. The mass accretion rate at inner disk can be derived using energy 

477 conservation law and Virial theorem. Assuming the disk radiation to be blackbody and considering 

478 both side of the disk, we get : 



m 



87rRl^km<^Tl,kev/^GMbh M^yr-^ (1) 



where a is the Boltzmann constant, G is the gravitational constant and T^in,keV is the inner 
disk temperature in keV and Rjn.fem is given by the equation 

Rin,km = V Ndiskbb X -Dio,fcpc/v^COSi (2) 

481 where ^diskbb is the normalization corresponds to the diskbb model component in XSPEC, Dio, fcpc 

482 is the distance to the source in the unit of 10 kpc and i is the disk inclination angle with respect 

483 to observer's line of sight. Inserting the equation (2) in (1), we get 



m 



'"''"''^ X iV^dtskbb X Dio^kpc/Vcosif (3) 



3GM 



The spectral analysis of a stable state (say, (j) class) in GRS 1915+105 using RXTE/PCA 
provided the average disk temperature of 1.25 + 0.02 keV and average disk normalization value of 
197.34 + 34.42. From the HIMS spectral analysis of IGR J17091-3624, we found the value of N diskbb 
to be 13.87 + 1.35 and the disk temperature to be 1.24 + 0.03 (See Table 2). Using this and assuming 
disk inclination angle (i) to be 70°, we found mGiJ5i9i5+i05/rii/Gi?,Ji709i-3624 ~ 3.78 + 0.15. If we 
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consider low disk inclination angle, then for i= 30", we found 1^1(3^51915+105/11170^717091-3624 ~ 
18.77 lb 1.23. Thus mass accretion rate is always higher in GRS 1915+105 compared to IGR 
J17091-3624, irrespective of the value of disk inclination angle. This is consistent with high X-ray 
flux observed in GRS 1915+105 than that seen in IGR J17091-3624 for same distance. We also 
calculated the relative viscous time-scale of the accretion flow using the following relation: 

= 4.3 X lO-^a-^m^'M-'Rl^,^ s (4) 

484 where t^^^ is the dynamic viscous time-scale, a is the viscosity parameter. We considered the value 

485 of a to be 0.01 for both sources since this is typical value for most of the BHXBs. Using relation 

486 ( 4) and above parameters, we found the value of ^'^is\GRSi9i5+W5/^vis\lGRJn09i-3624 *° ^ 0.676 

487 ± 0.018 when the disk inclination angle of IGR J17091-3624 is ~70". If the disk inclination angle 

488 is 30°, we found t'^^islGRSl9l5+l05/t^s\IGRJl709l-362A ^6 ~ 1.281 ± 0.177. Thus for a range of 

489 inclination angle between 30° and 70°, we found viscous timescale relative to GRS 1915+105 varied 

490 between 0.676—1.281. If we assume that 'p'-likc variability occurred in both sources at same viscous 

491 timescale then the disk inclination angle should be very small in IGR J17091-3624. 

492 Thus long term and continuous observation of this source in multi-wavelength band, particu- 

493 larly in hard X-rays as well as soft X-rays would reveal a lot of information regarding the radiation 

494 mechanism and accretion flow properties of BHXBs, and help to construct the complete picture 

495 of accretion flow around black hole systems by connecting GRS 1915+105 with other canonical 

496 BHXB systems. 



497 5. Conclusions 

498 Prom our analysis of SWIFT/XRT and RXTE/PCA data during 2011 outburst in IGR J17091- 

499 3624, and using the results from comparative study of timing and spectral properties with GRS 

500 1915+105, we conclude following points. 

501 1. We found that IGR J17091-3624 is the only LMXB which showed regular, repetitive and 

502 structured variability in the intensity pattern similar to the p-typc variability seen in GRS 

503 1915+105. Various parameters like burst profile, rise and decay profile, peak-to-minimum 

504 flux ratio, harmonics of mHz QPOs, 8—10 Hz QPOs and PDS show remarkable similarity. 

505 2. Several evidences like observation of variability, HID, PDS, power-law dominating spectra 

506 etc. provide some indications that the central compact object in IGR J17091-3624 may be a 

507 black hole. The detection of type-B and type-C QPOs along with p class activiity (as seen in 

508 GRS 1915+105 which is a black hole) is consistent with the black hole identification of the 

509 source. 

510 3. The source showed increase in flux while transiting from hard state to the HIMS. Later a 

511 transition from a typical irregular variability state, seen commonly in most of the BHXBs 
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512 to the regular, repetitive large variability state, seen only in GRS 1915+105 ('p'-type vari- 

513 ability). However, X-ray bursts occur in a state which is spectrally harder than that found 

514 in GRS 1915+105. Central black hole mass, its distance, absorption column density etc. 

515 may be important to explain the observed hardness ratio but at present we do not know 

516 exact answer. It is a topic for future investigations if these parameters or some unknown 

517 parameter /mechanism is responsible for higher hardness ratio during p class activity in ICR 

518 J17091-3624. 

519 4. Prom Table 1 and bottom panel of Fig. 9, it is clear that we detected several QPOs ranging 

520 from few mHz to 10 Hz or so. Such a large range in variability time-scale observed in 

521 this source indicate that the accretion disk which may be coupled with the corona, evolved 

522 significantly in this source. 

523 5. ICR J17091-3624 is the only source after GRS 1915+105 which show activities like '/9' class. 

524 In both sources, class activity is observed after the hard intermediate state where disk 

525 emission in the energy spectra is visible and contribute significantly to the total flux. GRS 

526 1915+105 never show typical low hard state where disk emission is not visible but IGR J17091- 

527 3624 shows low hard state where disk emission is not significant. This result agree well with 

528 the low hard state observations of other black hole X-ray binaries. Hence IGR J17091-3624, 

529 at certain spectral domain, behave similar to GRS 1915+105 whereas at some other spectral 

530 domain, it behaves like other canonical BHXBs. Hence, this source behaves as a missing link 

531 between GRS 1915+105 and other canonical BHXBs in understanding their observational 

532 dissimilarities, as well as accretion and radiation mechanisms. 
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Table 1: Observational details of IGR J17091-3624 using RXTE/PCA and/or SWIFT/XRT data 
during 2011 outburst. HIMS & TVS stand for Hard intermediate state and Intermediate variable 
state respectively 



ObHCrvatioii 

date (DD-MM-YY) 


MJD 


Iiistriiniont 


Observation 
ID 


Good time 
interval (sec) 


Avera}j;c 
soft color 


Average 
hard color 


lightcurve type and PDS features 


03-02-11 


55595.90 


XRT{PC) 


00031920001 


983.6 






iion-vaiiable, low flux and noise dominated (hard) 


09-02-11 




15601.06 


XRT(WT) 


00031921005 


1467.4 






low variability, low flux and noise dominated (hard) 


10-02-11 




55602. .58 


XRT(WT) 


00031921006 


2200.4 


- 


- 


low variability, low flux and noise dominated (hard) 


12-02-11 




55(10121 


XRTiWT) 


00031921008 


219S.1 






km' \";irial)ilit\", low flux and noisi* iloniinatcil fliardl 


lrl-li'^-11 




5.5()ll.5,21 


XRTi \\"T i 


l)003i!)21(}(}!l 


■''i-51.2 








14-02-11 




55606.16 


XRT(WT) 


00031921010 


1877.6 


- 


- 


low variability, low flux and noise dominated (hard) 


15-02-11 




55607.23 


XRT(WT) 


00031921011 


2101.6 






low variability, low flux and noise dominated (hard) 


16-02-11 




55608.23 


jYXv J. I vV L I 


q91 ni 9 
uuuo-i y ^ J- u J- ^ 


2182 3 






low variabihty, low flux and noise dominated (hard) 


18-02-11 




55610.17 


YRT^WT'l 
jYr\ J. I V V J. 1 


000*^1 Q91 ni ^ 


2176 5 






low variability, low flux and noise dominated (hard) 


20-02-11 




55612.25 


XRT(WT) 


00031921014 


2104.2 






low variability, low flux and noise dominated (haid) 


22-02-11 




55614.19 


XRT(WT) 


00031921015 


2050.2 






low variability, increased flux and noise dominated (hard) 


23-02-11 




55615.99 


PCA 


96065-03-01-03 


3184.4 


731 -1-0 01 3 


n Qfiq-l-n 021 


low vai'iahle high finx Hiid 4 29 IIz 
strong QPO foat\nc with weak red noise noticed (IIIMS) 


21-02-11 


55616.25 


XRT(WT) 


00031921016 


1065.6 






high flux and medium variability (IIIMS) 


26-02-11 


55619.01 


PCA 


QfiOfi '1-03-09-00 


1712.2 


U. 1 Ut:— 1— V7. UL/U 


n 400+0 01 Q 


strong QPO featine with weak white noise noticed (HIMS) 


02-03-11 


55622.59 
55622.52 


PCA 

XRT(WT) 


96420-01-01-00 
00031 991 OIK 

\.}\.}\.} ij ± 'J ^ 1-\J LO 


3744.3 
2340.6 


0.724±0.037 


0.390±0.078 


^ low variability, high fiiix and 4.18 IIz 

strong CJPO feature with weak w^hite noise noticed (HIIVIS) 


03-03-11 


55623.28 
55623.86 
55623.29 


PCA 
PCA 

XRT(WT) 


96420-01-01-01 
96420-01-01-02 
00031921019 


2064.5 
6944,4 
971.4 


0.684±0.034 
659+0 026 


0.364±0.023 
385+0 026 


~ 4.41 Hz & 9.07 Hz QPO feature and 

4 65 Hz strong QPO feature w^ith weak w^liite 
noise noticed; high irregular variability & high flux (HIMS) 


04-03-11 


55624.56 
55624.37 


PCA 

XRT(WT) 


96420-01-02-00 
00031991020 


1584.2 
2289.4 


0.623±0.028 


0.327±0.007 


^ 5.38 Hz strong QPO feature with 
low flux and noise dominated (HIM!S) 


06-03-11 


55626.39 


PCA 


96420-01-02-01 


65.5 


0.665±0.033 


0.349±0.014 


~ 5.53 Hz QPO feature with 

weak red noise noticed; low variability, 

low flux and noise dominated flllMSl 


l)7-li:;-ll 


.').".()_>7,ii2 


XliTi \\ T i 


l)00:li!)21u21 


2075.7 






luw \'arial )iUl \'. liij^h Ihix iuid ni-tise dnniiuiUrd (llhMSj 


08-03-11 


55628.47 
55628.23 


PCA 

XRT(WT) 


96420-01-02-02 
00031921022 


3202.2 
1704.1 


0.643±0.014 


0.356±0.005 


^ 5.39 Hz QPO feature with 

weak red noise noticed; low variability, 

low flux and noise dominated (HIMS) 


10-03-11 


55630.84 
55630.16 


PCA 

XRT(WT) 


96420-01-02-03 
00031921023 


1648.8 
1401.5 


0.669±0.019 


0.371±0.038 


^ 5.17 Hz strong QPO feature with 
weak red noise noticed; low variability, 
low flux and noise dominated (HIMS) 


12-03-11 


55632.54 
55632.73 


PCA 

XRT(WT) 


96420-01-03-00 

00031921024 


1520.4 

2182.3 


0.626±0.023 


0.328±0.024 


7.06 Hz QPO feature with weak red noise iiotieed: 
irregular variability, high flux and noise doiniiiated (HIMS) 


13-03-11 


55633.24 


XRT(WT) 


00031921025 


2009.8 






weak red noise noticed: irregular variability, high 
flux and noise dominated (IIIMS) 


14-03-11 


55634.72 


PCA 


96420-01-03-01 


3344.2 


0.532±0.022 


0.282±0.056 


-12.4 nillz and -2.59 IIz QPO features with 
very weak red noise notiecxl: variability 
become regular and repetitive (IVS) 


15-03-11 


55635.65 


XRT(WT) 


00031921026 


1481.2 






regular variabihty and high flux noticed (IVS) 


18-03-11 


55638.54 


XRT(WT) 


00031921029 


2164.6 






— 14.6 niHz QPO feature noticed 

low variability &: noise dominated (IVS) 


19-03-11 


55639.71 
55639.67 


PCA 

XRT(WT) 


96420-01-04-00 
00031921028 


2752.3 
2631.5 


0.502±0.007 


0.212±0.016 


^ 30.7 mHz strong QPO feature with 
strong red noise noticed (p) 


20-03-11 


55640.47 


XRT(WT) 


00031921030 


2352.5 






~ 30.2 mHz strong QPO with 
red noise noticed [p] 



Note. — Continuation of Table 1 
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Observation 

date (DD-MM-YY) 


MJD 


Instrument 


Observation 
ID 


Good time 
interval (sec) 


Average 
soft color 


Average 
hard color 


lightcurve type and PDS features 


22-03-11 


55642.20 
55642.35 


PCA 

XRT(WT) 


96420-01-04-02 
00031921031 


2976.4 
2260.3 


0.539±0.005 


0.294±0.089 


^^26.6 mHz and ~51.8 mHz QPO feature with 
strong red noise noticed (p) 


23-03-11 


55643.78 
55643.25 


PCA 

XRT(WT) 


96420-01-04-01 
00031921033 


1152.4 
643.5 


0.551±0.112 


0.325±0.030 


~ 12.3 mHz strong QPO feature with 
weak red noise noticed (p) 


24-03-11 


55644.75 


PCA 


96420-01-04-03 


2752.8 


0.571±0.028 


0.334±0.043 


^ 14.4 mHz & 8.06 Hz strong QPO feature with strong peaked 
noise at 0.06 Hz noticed; quasi regular, 
low variability (IVS) 


25-03-11 


55645.86 


PCA 


96420-01-0.5-02 


3312.3 


0.576±0.013 


0.371±0.()()4 


Weak red noise and a lentaUve QPO feature 

at ^ 6.98 Hz noticed qua^i-rogular, low variability (IVS) 


26-03-11 


55646.89 


XRT(WT) 


00031921034 


2160.1 






15.1 niHz QPO feature is noticed (p) 


27-03-11 


55648.01 


PCA 


96420-01-05-00 


32990.4 


0.522±0.037 


0.272±0.027 


~ 24.5 mHz & 8.03 Hz strong QPO feature with 

weak red noise noticed (/>) 


29-03-11 


55649.06 


PCA 


96420-01-05-03 


2640.3 


0.546±0.009 


0.331±0.023 


~24.8 niHz & ~49.6 niHz QPO feature with 

weak red noise noticed (/)) 


30-03-11 


55650.98 
55650.74 


PCA 

XRT(WT) 


96420-()l-().5-01 
00031921036 


2304.5 
331.65 


0.502±0.009 


0.279±0.019 


Weak red noise noticed & QPO feature 

al ~ 24.3 niHz, 45.2 niHz k 70.4 mHz is found (XRT)(p) 


31-03-11 


55651.88 


PCA 


96420-01-05-04 


3968.1 


0.513±0.024 


0.345±0.043 


QPO features at 28.1 mHz. 54.5 niHz. 84.6 mHz 
and 9.04 Hz are found witli weak red noise (p) 


02-04-11 


55653.70 


PCA 


96420-01-06-00 


3488.7 


0.502±0.019 


0.298±0.008 


Twin peak QPO like feature at 30.7 mHz & 37.1 mHz 
and strong QPO at 68.7 mHz witli 
strong red noise noticed {/i) 


03-04-11 


55654.89 


PCA 


96420-01-06-01 


6112.3 


0.532±0.008 


0.312±0.037 


QPO features at ~37.2 niHz, 75.3 mHz and --8.06 Hz 
with strong red noise noticed {p) 


05-04-11 


55656.71 


PCA 


96420-01-06-02 


5136.3 


0.540±0.010 


0.275±0.036 


QPO feature at '^40.4 mHz and 
strong BLN feature noticed (p) 


10-04-11 


55661.75 


PCA 


D6308-17-19-06R 


3376.6 


0.576±0.009 


0.283±0.027 


QPO features at ^^46. 5 mHz, '^94.0 mHz and 

~8.05 Hz QPO features are noticed with strong red noise (p) 
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Tabic 2: Best fit parameters (with la error bars) obtained by fitting tlie hard state and hard 
intermediate state spectra of IGR J17019-3624 with the powerlaw and the diskbb+powerlaw model 
respectively. 



Instrument 


Observation 


Date 


state/ 


r 


kTi„ 


^diskhb 




FluX,,o„,„(aii; 




x7<iof 


used 


ID 


(DD-MM-YY) 


class 




(keV) 












XRT 


00031921005 


09-02-11 


Hard 


-^■""-0.04 






0.601 






657.98/632 


XRT 


00031921006 


10-02-11 


Hard 


1 7c:+ll.05 
-^■'■^-0.04 






0.452 






473.58/632 


XRT 


00031921008 


12-02-11 


Hard 


-^■"•'-0.03 






0.872 






728.57/632 


XRT 


00031921009 


13-02-11 


Hard 


-'■"•''-0.04 






0.982 






667.64/632 


XRT 


00031921010 


14-02-11 


Hard 


1 70+0. 04 
-'■'•'-0.0,'i 






0.705 






694.01/632 


XRT 


00031921011 


15-02-11 


Hard 


-^■"'-0.04 






0.951 






726.74/632 


XRT 


00031921012 


16-02-11 


Hard 


-'■"^-O.o;! 






1.104 






529.34/632 


XRT 


00031921013 


18-02-11 


Hard 


-^■""-0.04 






1.035 






486.28/632 


XRT 


00031921014 


20-02-11 


Hard 


1 70+0.04 
-'■'^-0.0.5 






0.980 






509.43/632 


XRT 


00031921015 


22-02-11 


Hard 


r, r,r,+0.07 

^■^^-0.0.5 






2.472 






495.40/632 


PCA 


96065-03-01-03 


23-02-11 


IIIMS 


r, r,n+0.03 

^■^^-0.03 


-'■^"-o.oe 




0.791 


0.685 


0.106 


35.12/42 


PCA 


96065-03-02-00 


26-02-11 


IIIMS 


90+0.04 
^■■^^-0.07 


'^■''•'-0.06 


ii-snlt'o 


0.675 


0.597 


0.078 


43.59/42 


PCA 


96420-01-01-00 


02-03-11 


IIIMS 


99+0.02 
^■^^-0.03 


'■^"-0.04 


2-69;;;:;!; 


0.470 


0.427 


0.043 


44.52/42 


PCA 


96420-01-01-01 


03-03-11 


IIIMS 


9 r,„+0.08 

^■•'"-0.03 




4.36;;;:* 


0.492 


0.402 


0.09 


34.42/42 


PCA 


96420-01-01-02 


03-03-11 


IIIMS 


9 .< 1+0.04 
^- '^-0.03 


0.82;;;:;^ 




0.350 


0.329 


0.021 


40.78/42 


PCA 


96420-01-02-00 


04-03-11 


IIIMS 


9 oq+0.0.') 


1 90+0.03 

-^■^'^-0.03 


12.30;;;g 


1.812 


1.480 


0.330 


31.89/42 


PCA 


96420-01-02-01 


06-03-11 


HIMS 


44+0.0.') 


1 ir,+0.16 
-^■^^-0.12 


16.30;i;i 


1.828 


1.503 


0.33 


33.58/42 


PCA 


96420-01-02-02 


08-03-11 


IIIMS 


9 a+0.04 
^■'^'•'-0.03 


1 94+0.02 
-'■^*-0.03 


13.S7tYfl 


1.765 


1.360 


0.405 


40.75/42 


PCA 


96420-01-02-03 


10-03-11 


IIIMS 


9 oq+0.0.') 


1 90+0.03 
-^■^"^-0.04 


10.97;;- 


1.714 


1.422 


0.29 


41.17/42 


PCA 


96420-01-03-00 


12-03-11 


HIMS 


9 47+0.03 


0.891-- 




0.305 


0.281 


0.024 


42.12/42 



Note. — r is the power-law index, ^powerlaw is the power-law normalization, kTi„ is temperature of the inner disk 
edge in units of keV and 'Ndiskbb is the disk blackbody model normalization. Fluxtotoi is the total flux, Fluxdisktb is 
the flux of the disk component and Flwcpoweriaw is the flux of the power-law component. All the fluxes are calculated 
in the energy range 2.0-60.0 keV for combined PCA and SWIFT/XRT spectra and 0.3-10.0 keV for SWIFT/XKT 
spectra only. Fluxes are given in a unit of 10~^ ergs/s/cm^. 
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Table 3: Summary of Comparative study of some characteristic parameters of IGR J17091-3624 
and GRS 1915+105. 



Characteristic 


IGR J17091-3624 


GRS 1915-1-105 


parameters 


(variable state) 


Cp' class) 


Oscillation timescale 


~ 20 s - 50 s 


~ 40 s - 110 s 


Change in timescale 


increases with time 


increases with time 




~ 1.32 mHz/day 


~ 0.36 mHz/day 


Peak to minimum 


~ 2.91 - 3.82 


~ 3.38 - 4.28 


flux ratio 






X-ray flux and 


strong anti-correlation found 


strong anti-correlation found 


hard color relation 






Variability pattern 


slow exponential rise (ti=trise ~ 14.67) 


slow exponential rise {ti=trise ~ 13.68) 


(fltting with f(t)=ae*'''') 


fast decay (tj=t/ajj ~ -3.62) 


fast decay (ti=t/aj( ~ -2.48) 


QPOs 


28.1 mHz, 54.5 mHz (on 31 Mar. 11) 


9.7 mHz & 19.2 mHz (on 26 May 97) 




37.2 mHz, 75.3 mHz (on 03 Apr. 11) 


12.1 mHz & 24.3 mHz (on 03 Jun. 97) 




46.5 mHz, 94.0 mHz (on 10 Apr. 11) etc. 


18.5 mHz, 37.9 mHz (on 22 Jun. 97) 
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Fig. 1. — Overall nature of 2011 outburst of IGR J17091-3624. Top panel: plot of background- 
subtracted RXTE/PCA average count rate between 2.0-60.0 keV with time, middle panel: simul- 
taneous plot of background-subtracted SWIFT/XJ\T average count rate between 0.3-10.0 keV with 
time, Bottom panel: fractional RMS values of RXTE/PCA count rate as observed during hard in- 
termediate state and variable state. Due to contamination, RXTE/PCA has no observation when 
the source was in hard state. Vertical lines separate the hard state, hard intermediate state and 
the variable state. Intermediate variable states are shown by stars. 
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Fig. 2.— 0.3-10.0 keV filtered image of IGR J17091-3624 taken by tlie SWIFT/XRT in tlie photon 
counting mode on 03 February 11. 
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Fig. 3. — Result from SWIFT/'KBT observations of IGR J17091-3624 during the transition from 
the hard state to the hard intermediate state. Top left panel: 0.3-10.0 keV background-subtracted 
hght curves of the hard state and the hard intermediate state. Top right panel: 0.3-10.0 keV rms 
normalized power density spectrum of the hard state and the hard intermediate state. Bottom left 
panel: 0.6-7.2 keV unfolded energy spectrum, fitted model component (powerlaw) and residuals of 
the hard state. Bottom right panel: 0.6-7.2 keV unfolded energy spectrum, fitted model components 
(diskbb and powerlaw) and residuals of the hard intermediate state. 
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Fig. 4. — Result from observations during transition from low-variability class (hard intermediate 
state) to high variability class (variable state) via transition class (transition state) of IGR J17091- 
3624. Top left panel: 2.0-60.0 keV RXTE/FCA light curves of three types of classes. Top right 
panel: 2.0-60.0 keV RXTE/PCA rms normalized power density spectrum of three types of classes. 
Bottom left panel: evolution of classes from hard intermediate state to variable state in the Hardness 
Intensity Diagram. Bottom right panel: evolution of classes from hard intermediate state to variable 
state in the color color Diagram. 
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Fig. 5. — Evolution of variability pattern during variable state in IGR J17091-3624. All 5 panels 
show 2.0-60.0 keV background-subtracted PCA light curve for five different observations from 19 
March 11 to 10 April 11 with Is bintime. 
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Fig. 6. — Comparative study between the variable state variability class of IGR J17091-3624 and 
the 'p' class of GRS 1915+105. The left panel figures correspond to 2.0-60.0 keV RXTE/FCA 
lightcurve of IGR J17091-3624 on 31 March 11 (top) and 10 April 11 (bottom) and the right panel 
figures corresponds to 2.0-60.0 keV RXTE/PCA lightcurve of GRS 1915+105 on 26 May 97 (top) 
and 22 June 97 (bottom). Decrease in variability timescale with time was noticed in both sources. 
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Fig. 7. — Comparative study between the variable state variability class of IGR J17091-3624 and 
the V class of GRS 1915+105. Top panels show the 2.0-60.0 keV RXTE/PCA lightcurve of IGR 
J17091-3624 (left) and GRS 1915+105 (right). Bottom panels show the plot of hardness ratio 
(defined as the ratio of RXTE/PCA count rate between 12.0-60.0 keV and 2.0-12.0 keV) with time 
for IGR J17091-3624 (left) & GRS 1915+105 (right). Strong anti correlation between X-ray flux 
and hardness ratio was found in both sources. 
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Fig. 8. — Comparative study between the variable state variability class of IGR J 1709 1-3624 and 
the '/?' class of GRS 1915+105. Left panel: combined rise profile (left) and decay profile (right) of 
several bursts in IGR J17091-3624. Start of all rise profiles are normalized to sec and start of all 
decay profiles are normalized to 15 sec. An exponential growth function (ae^/*) fitted the rise profile 
better than a straight line with the significance of 0.312. An exponential decay function (ae"^/*) 
fitted the decay profile better than a straight line with the significance of 1—1.99 x 10~^. Right 
panel: combined rise profile (left) and decay profile (right) of several bursts in GRS 1915+105. 
Start of all rise profiles are normalized to sec and start of all decay profiles are normalized to 70 
sec. An exponential growth function (ae^/*) fitted the rise profile better than a straight line with 
the significance of 1—1.94 x 10"^". An exponential decay function (ae~^'/*) fitted the decay profile 
better than a straight line with the significance of 1—1.31 x 10^'^'^. 
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Fig. 9. — All left panel figures correspond to IGR J 1709 1-3624 and all right panel figures correspond 
to GRS 1915+105. Top panels: 2.0-60.0 keV RXTE/PCA rms normalized power density spectra. 
Upper Middle panels: Typical 8-10 Hz QPO in the PDS. Lower Middle panels: Hardness Intensity 
Diagrams for both sources for the same definitions of hard color and intensity. Both panels show 
that both sources evolved in a similar way in the PDS as well as in HIDs as the variability timescale 
decreases. Bottom panel: Evolutions of the characteristic frequencies of the variable state of IGR 
J17091-3624 (empty circles) and the p class of GRS 1915+105 (filled circles). In case of IGR J17091- 
3624, the time zero in x-axis corresponds to MJD 55635.37, and in case of GRS 1915+105, the time 
zero in x-axis corresponds to MJD 50594.49. This panel displays that within the same timescale 
of observations, IGR J 1709 1-3624 shows faster evolution of the variability frequency compared to 
the 'p' class frequency of GRS 1915+105. 



